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Spatial quantum enhancement effects are studied under a unified framework. It is shown that the multiphoton 
absorption rate of photons with a quantum-enhanced lithographic resolution is reduced, not enhanced, contrary 
to popular belief. The use of adiabatic soliton expansion is proposed to beat the standard quantum limit on 
the optical beam displacement accuracy, as well as to engineer an arbitrary multiphoton interference pattern for 
quantum lithography. The proposed scheme provides a conceptually simple method that works for an arbitrary 
number of photons. 
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In many optical imaging applications, such as atomic force 
microscopy |1] and nanoparticle detection |2j, precise mea- 
surements of the displacement of an optical beam are required. 
It is hence important to know what the fundamental limit on 
the accuracy of such measurements is placed by the laws of 
physics, and how one can approach this limit in an experiment. 
It is now known that if an optical beam consists of A^ indepen- 
dent photons with wavelength A, then the minimum uncer- 
tainty in its spatial displacement is on the order oiX/\/N, the 
so-called standard quantum limit |2|] . The ultimate uncertainty 
permissable by quantum mechanics, however, is smaller than 
the standard quantum limit by another factor of VA^ I13| . An 
experiment that beats this standard quantum limit with non- 
classical multimode light has already been demonstrated [4^. 
On the other hand, in other optical imaging applications, such 
as hthography, microscopy, and data storage, detection of ex- 
tremely small features of an object is desired. The feature 
size of an optical intensity pattern cannot be smaller than A, 
due to the resolution limit ||5|]. Multiphoton absorption allows 
detection of smaller feature sizes, and the minimum feature 
size of multiphoton absorption using a classical coherent light 
source is on the order of X/y/N |al, which can be regarded 
as the standard quantum limit on the multiphoton absorption 
feature size. Nonclassical light sources allow one to do better, 
and the ultimate limit is smaller than the standard one by an- 
other factor of \/N \6,i\. A proof-of-concept experiment of 
this resolution enhancement has also been demonstrated |a| . 
In the time domain, very similar quantum Umits on the posi- 
tion accuracy of an optical pulse can be derived fy]. Given the 
striking similarities among the spatiotemporal quantum limits, 
one expects them to be closely related to each other, yet the 
formalisms used to described each of them are vastly different 
IX^i^, so a more general formalism applicable to all spa- 
tiotemporal domains would greatly facilitate our understand- 
ing towards the spatiotemporal quantum enhancement effects. 

In this Letter, we apply the temporal formalism used by 
Giovannetti et at. J9|] to the spatial domain, and show that 
the uncertainty in the beam displacement and the spot size of 
multiphoton absorption are in fact closely related. Using this 
newly derived result, we demonstrate how arbitrary multipho- 
ton interference patterns can arise from a continuous super- 
position of coincident-momentum states. We further present 



an unfortunate result, namely, that the multiphoton absorp- 
tion rate is reduced if the quantum lithography resolution is 
enhanced, contrary to popular belief |6]. Finally, we take ad- 
vantage of the general spatiotemporal framework to show that 
the idea of adiabatic soliton expansion, previously proposed 
to beat the temporal standard quantum limit 1 10], can also be 
used to beat the standard quantum limit on the beam displace- 
ment accuracy, as well as generate an arbitrary multiphoton 
interference pattern, for an arbitrary number of photons. The 
use of solitons is an attractive alternative to the more conven- 
tional use of second-order nonlinearity for quantum informa- 
tion processing, because the soliton effect bounds the photons 
together and allows a much longer interaction length for sig- 
nificant quantum correlations to develop among the photons. 
Consider A^ photons with the same frequency co and polar- 
ization that propagate in the x — z plane. A general wavefunc- 
tion that describes such photons is given by II IB 
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where \k\,....kN) is the momentum eigenstate, ^i, ...,feAr spec- 
ify the transverse wave vectors of the photons along the x axis, 
and (p{ki,...,kN) is defined as the multiphoton momentum 
probability amplitude. The longitudinal wave vectors along 
the z axis are all assumed to be positive. Due to the resolu- 
tion limit, (p{ki,...,kN) is band-limited, i.e., ^{ki,...,kN) =0 
for any \ki\ > 2n/X. One can then define the corresponding 
quantities in real space. 
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, — 0(^i,...,^Ar)exp(/feiXi + ... + /^ArJt:Ar), (3) 

2n \J2% 

1^) = —?= / dxi...dxt^\^{xx,...,xj^)\xx,...,xtj), (4) 

where i//(xi, ...,XAr) is the multiphoton spatial probability am- 
plitude. and y/ are subject to normalization conditions 



J dki...dkN\^\^ — J dxi...dxN\\i/\^ — 1, and (j) and Xj/ must be 
symmetric under any exchange of labels due to the bosonic 
nature of photons. The magnitude squared of i^ gives the joint 
probability distribution of the positions of the photons, 

{■.Iixi)...IixN) :) - ^(^|At(xi)...A^(x^)A(xi)...A(x^)|^) 

(5) 
(6) 



\\j/{xi,...,xn)\^, 



where A{xi) and A^(x,) are the spatial annihilation and cre- 
ation operators respectively. The statistical interpretation of 
y/ is valid because we only consider photons that propagate 
in the positive z direction. The above definition of a multi- 
photon state is more general than those used by other authors, 
in the sense that we allow photons with continuous momenta, 
compared with the use of only one even spatial mode and one 
odd mode by Fabre et al. (3|], the use of only two discrete mo- 
mentum states by Boto et al. [fl], and the use of many discrete 
momentum states by Bjork et al. |7|]. 

The displacement of an optical beam can be represented 
by the operator X = {I /N) J dx xA'' {x)A{x). Applying X to 
\xi,...,xn) gives 



X\x, 



,xn) 




,xn), 



(7) 



so the beam displacement is, intuitively, the mean position of 
the photons under the statistical interpretation. If we assume 
that (x) = for simplicity, the root-mean-square displace- 
ment uncertainty is given by 



-/^2\l/2 



AX = (X^) 
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It is often more convenient to use a different system of coor- 
dinates as follows 1I12I1 . 
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X,i^\,...,N-\, ^N = -Y.^i- 

i=l 

(10) 



X is therefore the "center-of-mass" coordinate that character- 
izes the overall displacement of the optical beam, and ^,'s are 
relative coordinates. Defining a new probability amplitude in 
terms of these coordinates. 



xi/'iX,^u-,^N-i) = ViX + ^u-,X + ^^ 



(11) 



we obtain the following expression for the displacement un- 
certainty. 



AX = 
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which is the marginal width of i// with respect to X. 

On the other hand, the dosing operator of A^-photon absorp- 
tion is given by 



(:/^W:)c<|v/(x,x,...,x)|2 = |v/'(x,0,...,0)|^ 



(13) 



which is, intuitively, the probability distribution of all A^ pho- 
tons arriving at the same place x. Hence, designing a speci- 
fied multiphoton interference pattern in quantum lithography 
is equivalent to engineering the conditional probability distri- 
bution |v/'(X,0,...,0)|2. 

In particular, the spot size of multiphoton absorption is the 
conditional width of Xj/' with respect to X, 
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Despite the subtle difference between the marginal width and 
the conditional width, if y/' can be made separable in the fol- 
lowing way. 



V/'(X,^i,...,^A,_i) = V/(X)v/,.w(<^i,...,<^A,_i), 



(16) 



then both widths are identical, and one can optimize the mul- 
tiphoton state simultaneously for both applications. 

The standard quantum limit on the uncertainty in X is ob- 
tained when the photons are spatially independent, such that 
V/(xi,...,XAr) =/(xi).../(x7v). For example, if /(x) is aGaus- 
sian given by f{x) °^ exp (^—K^x^/2), then both the marginal 
and conditional uncertainties in X are 
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Similar to the optimization of temporal position accuracy |9|], 
the ultimate quantum limits on spatial displacement accuracy 
and quantum lithography feature size are achieved with the 
following nonclassical state. 



^) = dkG{k)\k,k,...,k). 



(18) 



The momentum probability amplitude is then 



^iki,...,kN)^Giki)d{ki-k2)d{ki~k3)...d{ki-kN), 

(19) 

which characterizes A^ photons with coincident momentum. 
The spatial amplitude is thus given by 

dk 



xi/{X,^u-,^N-i)= I -^G{k)cxp{iNkX)=g{NX), 



(20) 



which is a function of X only and can be understood as a 
continuous superposition of A^-photon coincident-momentum 
eigenstates, each with an effective de Broglie wavelength 



equal to 2n/{Nk). This representation is equivalent to Boto et. 
fl/.'s proposal |6], when G{k) ~ 5{k~ ko) + 5{k + ko), where 
^0 is the transverse wave vector of either arm of the inter- 
ferometric scheme. The multiphoton interference pattern is 
therefore trivially given by \g{NX)\-^, the magnitude squared 
of the Fourier transform of G{k). An arbitrary interference 
pattern can hence be generated, if an appropriate G{k) can be 
engineered. This approach of designing the multiphoton in- 
terference pattern should be compared with the less direct ap- 
proaches by the use of discrete momentum states latzll- With 
the resolution limit, G{k) — for |A:| > In/X, so, given the 
Fourier transform relation between G{k) and g{NX), the min- 
imum feature size of multiphoton interference is on the order 
ofX/N. 

To compare the ultimate uncertainty in X with the stan- 
dard quantum limit, let G{k) be a Gaussian given by G{k) oc 
exp [— fc^/(2K"^)] , then the uncertainty in X becomes 
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(21) 



which is smaller than the standard quantum limit, Eq. Hit , by 
another factor of VA^, as expected. 

Let us recall Boto et al.'s heuristic argument concerning 
the multiphoton absorption rate of entangled photons. They 
claim that, because entangled photons tend to arrive at the 
same place at the same time, the multiphoton absorption rate 
must be enhanced |6]. If photons tend to arrive at the same 
place, then the uncertainties A(^,- in their relative positions ^, 
must be small. However, the spatial probability amplitude that 
achieves the ultimate lithographic resolution, Eq. ( I20> . is a 
function of X only, which means that the uncertainties in ^, 
are actually infinite. In general, any enhancement of resolu- 
tion with respect to X must result in increased uncertainties in 
the relative positions £,{, in order to maintain the same max- 
imum bandwidth. Hence, Boto et a/.'s argument, as far as 
the spatial domain is concerned, manifestly does not hold for 
photons with a quantum-enhanced lithographic resolution. In 
fact, the opposite is true: although these photons have a re- 
duced uncertainty in their average position X, they do not ar- 
rive at the same place as often, and the multiphoton absorption 
rate must be reduced. 

To observe this fact, consider the total multiphoton absorp- 
tion rate 



dx{:f{x):)oc dX\\i/'{X,Q,...,Q)\^ 



(22) 
Because \j/' must satisfy the normaUzation condition, 

^JdXd^i...d^N-i\w'{X,^u-AN-i)\^^h (23) 

JdX\\j/(X,0,...,0)\^ is inversely proportional to (A^,)'^"^ 
An increase in each A^, by a factor of y hence reduces the 
total absorption rate by a factor of y^^ ' . 

With all that said, if the multiphoton absorption rate is re- 
duced due to a quantum-enhanced lithographic resolution, one 



can still compensate for this rate reduction by reducing the rel- 
ative temporal position uncertainties of the photons 11311 . 

We now turn to the problem of producing the nonclassi- 
cal multiphoton states for spatial quantum enhancement by 
the use of optical solitons. Consider the Hamiltonian that de- 
scribes the one-dimensional diffraction effect and Kerr non- 
linearity on an optical beam in a planar waveguide. 
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(24) 



where b is the Fresnel diffraction coefficient, assumed to be 
positive, and c is the negative Kerr coefficient, assumed to be 
negative, so that b/c < and solitons can exist under the self- 
focusing effect. The soliton solution of the spatial amplitude 
for A^ photons under this Hamiltonian is 11411 
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(25) 



where C = ^J{N -ly^c/by^-^ /{2n) and G{k) is determined 
by the initial conditions. If initially the photons are uncorre- 
cted, G{k) can be approximated as lirH 



G{k) oc exp 
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(26) 



where ^ is a parameter on the order of unity JlOll . and Wo is the 
initial soliton beam width. The probability amplitude can be 
written in terms of the center-of-mass and relative coordinate 
system defined in Eqs. dlOt as 



V/' = C 



dk 



2k 

2, 



G(^)exp 



iNkX + -Y^\^,-^j\ 



ibNk'-t + i-^NiN^-~\)t 
I2b ^ ' 



(27) 



which is separable in the way described by Eq. ( I16> . meaning 
that the conditional width and marginal width with respect to 
X are identical. 

If we adiabatically reduce c or increase b along the waveg- 
uide by, for example, increasing the waveguide modal thick- 
ness, then we can reduce the uncertainty in the relative mo- 
menta of the photons and increase the uncertainty in the rel- 
ative positions |10]. Classically, we expect the soliton beam 
width to expand and the spatial bandwidth to be reduced. But 
the most crucial difference in the quantum picture is that the 
center-of-mass coordinate X remains unaffected during the 
adiabatic soliton expansion, apart from the quantum disper- 
sion term —ibNk^t. In the limit of vanishing c/b, the wave- 
function approaches the ultimate multiphoton state given by 
Eqs. (d, (El, and (|20j. 

As pointed out in Ref. |10], the quantum dispersion term 
can be compensated if the soliton propagates in a second 



medium with an opposite diffraction coefficient b' . Full com- 
pensation is realized when /q b{t)dt = — Jq b'{t)dt, where T 
is the propagation time in the first medium and T' is the prop- 
agation time in the second medium. Negative refraction 11511 
is hence required in the second medium. Ideally the second 
medium should also have a Kerr coefficient c' opposite to the 
final value of c in the first medium, such that c' = —c(r), so 
that b' /c' < and the quantum soliton maintains its shape, but 
in practice c' — also suffices, because the multiphoton spec- 
trum 1 01^ remains unchanged in a linear dispersive medium 
while the quantum dispersion is being compensated. 



effects, such as loss, become significant. However, for quan- 
tum lithography, an ideal high-number-photon absorption ma- 
terial is difficult to obtain, so a giant Kerr nonlinearity, such 
as that theorized in a coherent medium Il6ll . is required to 
produce a few-photon soliton. 

In conclusion, spatial quantum enhancement effects are 
studied under a general framework. It is shown that the mul- 
tiphoton absorption rate is reduced if the lithographic resolu- 
tion is enhanced. The use of adiabatic soliton expansion is 
proposed to beat the spatial standard quantum limits for an 
arbitrary number of photons. 
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FIG. 1: (color online). Schematics of the spatial quantum enhance- 
ment setup via adiabatic soliton expansion. 

So far we have worked in the paraxial regime, so a 4/ 
imaging system with spatial phase modulation in the Fourier 
plane can effectively mimic the behavior of negative refrac- 
tion. Consider the system depicted in Fig. [2 After the beam 
goes through the nonlinear medium and the first Fourier lens, 
in the Fourier plane, the multiphoton amplitude becomes 



V/(xi 



,xn) 













(28) 



A quadratic spatial phase modulation in the Fourier plane, by 
a lens for example, can hence act as negative refraction in the 
paraxial regime and cancel the quantum dispersion term. 

In general, Fourier-domain modulation with a transfer func- 
tion H{k) can be used to shape G{k), if the ultimate multi- 
photon state is achieved, because all photons have coincident 
momenta, resulting in an output given by G{k)H^ (k). A de- 
sired multiphoton interference pattern \g{NX)\-^ can hence be 
engineered by spatial modulation in the Fourier plane. 

If the ultimate multiphoton state is achieved via adiabatic 
soUton expansion, the bandwidth of the optical beam is the 
same as that of G{k), which is on the order of l/( v^Wo), and 
is much smaller than the resolution limit. The second Fourier 
lens in Fig. [l] should therefore have a small focal length /2 
to demagnify the optical beam and increase the bandwidth of 
G{k). 

Current technology should be able to expand a spatial soli- 
ton with A^ ^ lO'" photons by a few times before decoherence 
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